The effects of time delays in adaptive phase measurements 
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It is not possible to make measurements of the phase of an optical mode using linear optics without 
introducing an extra phase uncertainty. This extra phase variance is quite large for heterodyne 
measurements, however it is possible to reduce it to the theoretical limit of log n/(4n 2 ) using adaptive 
measurements. These measurements are quite sensitive to experimental inaccuracies, especially time 
delays and inefficient detectors. Here it is shown that the minimum introduced phase variance when 
there is a time delay of r is r/(8n). This result is verified numerically, showing that the phase 
variance introduced approaches this limit for most of the adaptive schemes using the best final 
phase estimate. The main exception is the adaptive mark II scheme with simplified feedback, which 
is extremely sensitive to time delays. The extra phase variance due to time delays is considered for 
the mark I case with simplified feedback, verifying the t/2 result obtained by Wiseman and Killip 
both by a more rigorous analytic technique and numerically. 
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It is well known that any measurement of the phase of an electromagnetic field using linear optics will have 
, an uncertainty that is greater than the intrinsic quantum phase uncertainty of the state Q. The standard phase 
measurement method is the heterodyne scheme, where the signal is combined with a strong local oscillator field with 
^\ , a slightly different frequency. This means that all quadratures of the field are sampled approximately equally. This 
phase measurement method introduces a phase variance of approximately I /(An), where n is the mean photon number 
O ; of the field. 

This is not so significant for measurements on coherent states, where it is the same size as the intrinsic phase 
uncertainty of the state. For states with reduced phase uncertainty, however, the uncertainty in the measurement 
will be far greater than the intrinsic phase uncertainty. For example, minimum phase uncertainty states have a phase 
variance that scales as l/(n 2 ) [^J. 

It is possible to dramatically improve on heterodyne measurements by using a local oscillator phase of$ = <y3 + 7r/2, 
where ip is the signal phase (homodyne detection). This has the drawback that the phase must be known in advance. 
Adaptive phase measurements (l], |3|, |^, ^| attempt to approximate a homodyne phase measurement by adjusting 
the local oscillator phase based on data obtained during the measurement. 

There are several different variations of adaptive phase measurements. The adaptive mark I scheme gives improved 
^J^" results for small photon numbers, but worse results for large photon numbers. The adaptive mark II scheme gives an 
introduced phase variance of l/(8n 15 ) [^, a significant improvement over heterodyne measurements. By a more 
• i-H \ sophisticated feedback algorithm it is even possible to obtain the theoretical limit of logn/(4n 2 ) 

These adaptive measurement schemes are sensitive to experimental imperfections, most notably imperfect detectors 
and time delays. The effect of imperfect detectors is fairly straightforward, introducing a phase variance of approxi- 
mately (1 — ?j)/(4rjn) j|. The effect of time delays is more difficult to estimate. Some highly simplified calculations 
indicate that the excess phase variance due to time delays is t/2 for mark I measurements (where r is the time delay), 
and r/(2n) for mark II measurements 0. 

Here we repeat these derivations more accurately, and show that while the result for mark I measurements is 
reasonably accurate, the perturbation approach is inadequate to obtain a consistent result for mark II measurements. 
We consider an alternative derivation that gives the minimum phase variance when there is a time delay. In section ^ 
we evaluate the phase variance with time delays numerically and show that for most of the measurement schemes the 
phase variance approaches this limit for large time delays. 

2. BACKGROUND THEORY 

Before we proceed to determining the effect of time delays, we will briefly outline the background theory for adaptive 
measurements. For more details see jl], Q. First the photocurrent for dyne detection Jlj is given by 

I(v)dv = 2Re(a 1 ,e- J<E, ^)dv + dW(v). (1) 
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Here a v is the scaled coherent amplitude of the signal, v is the scaled time and is the local oscillator phase. The 
systematic variation with time of the coherent amplitude due to the mode shape is scaled out, and time is scaled to 
the unit interval. We also define the variables 



A v = / 7(w)e l * (u) du, (2) 
Jo 

B v =- f e 2 ^du, (3) 
Jo 

C v = A v v + B V A* V . (4) 



We omit the subscripts to indicate final values. The best estimate of the phase at time v is given by argC„, and if 
\B V \ is small then aigA v is also a good phase estimate. 

For mark I measurements aigA v is used as the intermediate phase estimate (so the local oscillator phase is $ = 
arg A v + 7r/2) and also as the phase estimate at the end of the measurement. For mark II measurements we use the 
same intermediate phase estimate but we use the best estimate argC at the end of the measurement. In it is 
shown that using arg A v as the intermediate phase estimate is equivalent to varying the local oscillator phase as 

k I(v)dv ,_. 

Using the feedback in this form allows us to use a simplified feedback circuit experimentally. This feedback is no 
longer equivalent to using a phase estimate of &vgA v when there are time delays in the system, and we therefore 
consider the cases with arg A v feedback and simplified feedback separately. 

The next feedback scheme uses a phase estimate intermediate between arg A v and arg C v , specifically 

V = Mg{A%C l v -*). (6) 

The simplest case is where e does not depend on time. For each mean photon number there is an optimum value to 
use, and provided these optimum values are used this method gives better results than mark II measurements || . 
We can get very close to the theoretical limit of logn/(4n 2 ) if we use a time dependent e given by 



v*- b v * 

Briefly explaining the reason for the theoretical limit, the probability distribution for A and B is proportional to 
\{(3, CIV')! 2 ! where \ip) is the signal state and |/3, C) is a squeezed state 

|/3,C) = e X p(/3at-/3*a)exp[iC*a 2 -K( at ) 2 ]l°)» ( 8 ) 

where 

C_ 
l~B 2 



^ = 7-^ (9) 
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This means that the phase variance introduced is approximately equal to the phase variance of the squeezed state 
Q. The photon number of this squeezed state will be approximately the same as that of the input state, so the 

theoretical minimum introduced phase variance is that of an optimized squeezed state with photon number n. As 

shown in 0, this scales as logn/n 2 . 

The phase measurement scheme with time dependent e (0) is not quite at the theoretical limit because it produces 

squeezed states \/3, £) that are slightly more highly squeezed than optimum. In |g) we show how this can be corrected, 

but we will not consider the case with these corrections here because the time delays cause the squeezed state to be 

less squeezed than optimum anyway, so these corrections are not needed. 



3. PERTURBATION APPROACH 



Mark I 



Now we estimate the effect of time delays on simplified mark I measurements in a similar way to that done in [Q , 
but using fewer of the simplifications used there. Without a time delay the stochastic differential equation (SDE) for 



the phase estimate is 



I(v)dv 



1/2 [-2a sin tp v dv + dW(i 



(11) 



In this expression we have taken the input phase to be zero. For some time v\ the phase will come to lie near 0, so 
we linearise around tp v = 0. The result, which will be valid for v\ < v < 1 is 



d(p v = v~ 1/2 [-2aip v dv + dW{v)]. 

Including the time delay the SDE is 

dtp v = v~ 1/2 [-2aip v ^ T dv + dW(v)]. 
Now we treat the delay perturbatively. We write the solution to the perturbed equation as 

<p v =tpW +otT 0W +0(a 2 T 2 ). 
The zeroth-order term obeys the SDE for no delay (|l2|), so the first-order correction obeys 

ard^ 1 ' = 2av~ 1 ' 2 {ip^ - $± T )dv - 2a 2 rv- 1 / 2 ^ v 1) _ T dv. 

Therefore to first order in r we have 

d^ 1 ) = 2v- 1 / 2 dip^ - 2av- 1 / 2 ^dv 

= 2v- 1 / 2 {v- 1 ' 2 [-2a^ ) dv + dW(v)}} - 2av~ 1/2 ip { v 1) dv 

= -2av~ 1/2 ^dv ~ — ip^hlv + -dW(v). 

V V 

It is straightforward to show that the solution to the zeroth-order equation is 



(0) _ iaiV^-WfiiO) , / e 



Using this in equation ([l6]) and multiplying on both sides by e Aa ^" gives 

4a 



d(e 4Q ^^ 1} ) = 
Integrating then gives the solution 



e 4oV5T^(0) + 



V 4ay/u 



dW(u) 



dW{u). 



dv + -e 4a ^dW(v). 



rv c rvi 



v 4a 



s e 4a( v / M-\A') 



dW(u) 



4a(v^T-v^)^(0) 



ds+ [ -e 4a( ^-^dW(s). 
J Vl s 



(12) 
(13) 
(14) 

(15) 



(16) 



(17) 



(18) 



(19) 



In this approximation the mark I phase estimate is given by ip\ = tp^ + anp^ . To first order in r and a 1 this has 
a variance of 



1 



(20) 



Here we have used the known variance of l/(4a) of the zeroth-order term. Evaluating the second term on the 
right-hand side we find 



((^W) - e- 8 ^ 1 -^) (<^«) -Malog-Me-M 1 ^) ^(0)2^ 



4a f s e 8a (Vu-i) 



duds 



1 2e 8a(v^-l) 



,1.5 



-du. 



(21) 
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The first two terms decrease exponentially with a and may therefore be omitted. Exchanging the order of the integrals 
in the third term and integrating gives 



-(o) -(i) 



4a 



1 logu 



1 2 e 8a<y«-l) 



,1.5 



du. 



Now we change variables to s = 1 — ^/u, so du = —2(1 — s)ds. Then we obtain 



£?M ly ) = 16a 



Expanding in a Maclaurin series in s gives 



-(o) -(1) 



x -^ log(l- a ) Sai 

(1-8) 



ds + 4 



1~VW g-8as 



ds. 



16a 

l- 



(i 



-16a 



+ 4 



1 

8a + 



^0 + \s 2 + 0(s 3 ))e- 8as ds 
1 (1 + 2s + 3s 2 + 0(s 3 ))e~ 8as ds 
' 4 +0(a" 4 ) 

1 0(a- 4 



ia) 2 (8a) 3 
2 6 



4a 



a) 2 (8a) 3 
0(a- 2 ). 



(22) 



(23) 



(24) 



Note that the upper bound at 1 — ^/vi has no effect since it gives a term that decays exponentially with a. Thus we 
find that the total phase variance is 



(25) 



This provides a good verification of the result obtained by the highly simplified method in |J. 

Note that this result is based on continuing to use the intermediate phase estimate at the end of the measurement. 
If we use the phase estimate &rgA at the end of the measurement we will get a different result, one that we cannot 
predict using this approach. 



B. Mark II 



If we try to use the same approach for the mark II case it does not seem to be possible to obtain a consistent result. 
To illustrate this, we will briefly outline the derivation. From Q, the mark II phase estimate is effectively a time 
average of the mark I phase estimates: 

l 

tftdt. (26) 

In order for this to be consistent with the above theory we will take the average only from time v\, then take the 
limit of small v\. In perturbation theory the mark II phase estimate is 



$°>+or#]di. (27) 



We find that the variance is 



(4) = f dt f df U ^^) + 2ar f dt f dt> U a) ^ } ) ■ (28) 



The first term is fairly well behaved, and in the limit of small v\ we get 



Cdtfdt'^^)^±,. (29) 
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Simplifying this, the third and fourth terms cancel, and the first two terms give 



Unlike the results for the mark I case, all the terms depend on the conditions at time v\. This means that it is not 
possible to obtain an unambiguous result in the limit v\ — » 0. Therefore we consider an alternative approach for 
estimating the increase in the phase variance due to the time delay. 



4. THEORETICAL MINIMUM 



The alternative method of obtaining an estimate of the time delay is to consider the squeezed state |/3, Q in 
the probability distribution. As was explained above, the excess phase variance due to the measurement scheme is 
approximately the phase variance of this squeezed state. 

From 0, the phase variance of a squeezed state is given by 

(A0 2 )«^±i + 2crfc(y2^), (32) 

where no = n p e 2< » for real £. Here we use the subscript p to indicate the mean photon number of the squeezed state 
in the probability distribution, as opposed to the photon number of the input state. The average value of n p will be 
close to the photon number of the input state. 

For states that are significantly less squeezed than optimum, the second term is negligible and we can omit the 
term of order n~ 2 . Then this simplifies to 

<A^£. (33, 

Since n p will be close to the photon number of the input state, it is reasonable to replace it with n. 

When we have a delay of r in the system, before time r we have no information about the phase of the system 
to use to adjust the local oscillator phase. Therefore we must use a heterodyne scheme for this time period, rapidly 
varying the local oscillator phase. This means that B T will be equal to zero, and no matter how good the phase 
estimate is after time r, the largest the magnitude of B v can be made is v — t. Then at the end of the measurement, 
the largest \B\ can be is 1 — t, and the largest |£| can be is atanh(l — r). 

The lower limit to the phase variance of £) when there is a time delay of r is therefore 

„ — 2atanh(l — r) 

<^ 2 > min Tn ■ (34) 

pa (35) 

This is therefore also the lower limit to the introduced phase variance when there is a time delay of r. We can expect 
that the introduced phase variance will be close to this for states of small intrinsic phase variance, as there will quickly 
be very good phase estimates available for the feedback phase. In addition the time delay must be sufficiently large 
that the phase variance given by this expression is significantly above the introduced phase variance for no time delay. 

This result obeys the same scaling law as the result given in EQ, but it is a factor of four times smaller. Note, 
however, that the limit condition for the result given in is that to. is small, whereas the above result should only 
be accurate when both a and r are reasonably large. The result here also differs in that it is the limit for the total 
introduced phase variance, rather than just the extra phase variance due to the time delay. 
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5. NUMERICAL RESULTS 

These analytic results were also tested numerically. The numerical techniques used were similar to those used 
in reference rol. Minimum uncertainty squeezed states were used, with the stochastic differential equations for the 
squeezing parameters H as given in 0. For all calculations 2 20 time steps were used, and calculations were performed 
with time delays of 2™ time steps, where n varies from to 18. 

For the first 2 n time steps the local oscillator phase was rotated by ir/2 each step. For the following time steps the 
data up to the time step 2™ before the current time step was used. For a delay of 2° = 1 time steps the data from 
the previous step is used, corresponding to the technique for no time delay. 

Numerical results for four different phase feedback schemes were obtained: 

(a) The simplified feedback for mark I and II measurements, where 

dp v = =— . (36) 

(b) The unsimplified feedback, where the phase estimate is 

tp(v)=axgA v . (37) 

(c) The phase estimate that is intermediate between arg A v and the best phase estimate 

0(«) = arg(^Cj-«), (38) 

where e is a constant. 

(d) The same as in (c), except that the value of e varies with time as 



v 2 -\b v \ 2 r~r~ 
e(v) = T^Vl^- (39) 

A. Comparison with perturbative theory 

First we consider the case of simplified feedback, and consider th e var iance in the final value of the feedback phase, 



rather than the phase of A or C. This case was examined in section 3 A , and the extra phase variance due to the time 
delay is r/2 according to that analysis. The extra phase variance is plotted for four different mean photon numbers 
in figure \U. 

In determining the extra phase variance due to the time delay, an estimate must be made of the phase variance 
with no time delay. For the results shown in figure R], the phase variances are very close for the first six or so time 
delays. The estimate of the phase variance with no time delay was taken to be the minimum of these results. 

The theoretical value of r/2 is also plotted in figure [l]. As can be seen, many of the results are close to the theoretical 
line for the intermediate time delays. For small time delays the extra phase variance due to the time delay is too 
small a fraction of the total phase variance for the results to be accurate. The reason why the results deviate from 
the theoretical result for large time delays is that this theoretical result is for the limit of small ar. Note also that the 
results for larger photon numbers deviate from the theoretical result for smaller r than the results for smaller photon 
numbers. This is also to be expected from this limit condition. 

We also used fitting techniques to determine how closely the numerical results agree with the theoretical value. 
A linear fit of the phase variances against r was performed, and data for ar above about 0.3 was omitted, as this 
was where the results started to increase dramatically. The average slope obtained was 0.39 ± 0.06, in reasonable 
agreement with the theoretical value of 0.5. 

As was mentioned above, the result for the additional phase variance due to the time delay is only valid for the 
variance in the final value of the phase estimate, which is not the same as arg A when there is a time delay. In figure |^ 
we have plotted the variation of the phase variance with time delay for three alternative final phase estimates, tp, 
arg A and argC. This is for a photon number of approximately 332 000, and is fairly representative of the results for 
other photon numbers. 

As can be seen, for very small time delays the variances in the (p and arg A phase estimates are almost identical. As 
the time delay is increased, however, the variance of ip increases, but the variance of arg A decreases. This is because, 
as the intermediate phase estimate gets worse, the value of \B\ decreases. This means that A is closer to C, so arg A 
is closer to the best phase estimate. Note, however, that the variance of arg A rises again, and does not converge to 
argC for large time delays. This is because \B\ does not fall to zero. 
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FIG. 1: The extra phase variance (in the final value of the intermediate phase estimate for the simplified feedback) due to the 
time delay plotted as a function of time delay for four different mean photon numbers. The data for a mean photon number of 
121.590 is shown as crosses, for a photon number of 1 576.55 as circles, for a photon number of 22 254.8 as asterisks and for a 
photon number of 332 067 as plusses. The theoretical value of r/2 is plotted as the continuous line. 
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FIG. 2: The phase variance of three alternative final phase estimates for simplified feedback with a time delay plotted as a 
function of time delay. The results for the final value of the intermediate phase estimate are plotted as a continuous line, for 
arg A as a dotted line, and for arg C as a dash-dot line. All results are for a photon number of 332 067. 
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FIG. 3: The phase variance introduced for three different phase feedback schemes plotted as a function of time delay. The 
dotted line is for simplified feedback, the dash-dot line is for the corrected simplified feedback, and the circles are for unsimplified 
arg A v feedback. The best phase estimate arg C is used in all three cases. The continuous horizontal line is the phase variance 
for heterodyne measurements, and the continuous diagonal line is the theoretical limit. All results are for a photon number of 
332 067. 



B. Comparison with theoretical minimum 



Now we consider the variance in the phase of C. As was explained above, the theoretical lower limit to the 
introduced phase variance is r/(8n). We have plotted the introduced variance in the best phase estimate argC and 
the theoretical limit in figure pi For additional accuracy we have plotted 



-2atanh(l-r) 

An 



(40) 



as this will continue to be accurate for time delays that are a large fraction of 1. This plot is for a photon number 
of 332 000, and similar results are obtained for other photon numbers. As can be seen, the phase variance is well 
above the theoretical limit. For large time delays the phase variance approximately converges to the heterodyne phase 
variance, also shown in figure |^. 

The phase variance introduced for mark II measurements with the unsimplified aigA v feedback is also shown in 
figure |^. The phase variance introduced for this case increases far more slowly with the time delay, and for larger 
time delays it is very close to the theoretical limit. These results indicate that if there is any significant time delay in 
the system the simplified feedback will give a far worse result than using arg A v . 

It is possible to make a correction to the simplified phase feedback scheme that improves this result somewhat. 
Many different alternatives were tried, and the one that gave the best results was 



d0 v = 



I(v)dv 
y/v + ar 



(41) 



This correction is based on the fact that \A V \ is larger than y/v when the phase estimate is worse than aigA v . (From 
PL the factor of y/v in the simplified feedback comes from a factor of \A V \.) 

The results for this correction are also shown in figure |[ The phase variances obtained in this case are significantly 
below those for the plain simplified feedback, but are still far above the results for the unsimplified aigA v feedback. 

Now we consider the results for better intermediate phase estimates that are between arg A v and arg C v . The phase 
variance introduced for the constant e case and the theoretical limit are shown in figure These results are again 
for a photon number of about 332 000. The results for this case are even closer to the theoretical limit than those for 
the mark II case. 
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FIG. 4: The phase variance introduced for better intermediate phase estimates plotted as a function of time delay. The plusses 
are for the constant e case and the crosses are for the time-dependent e case. The theoretical limit estimated using the mean 
inverse photon numbers obtained from the time-dependent e case is plotted as the dotted line, and the theoretical limit using 
the input photon number is shown as the continuous line. All results are for a photon number of 332 067. 

The phase variance introduced for the feedback with time dependent e is also plotted in figure ||. The results for 
this case converge to the theoretical limit at smaller time delays than for the constant e case. For the larger time 
delays the results for the two cases are about the same, slightly above the theoretical limit. 

In both cases we find that for large time delays the phase variance is still noticeably above the theoretical limit, and 
that the values of B obtained are too close to 1 — r to account for this difference. This difference appears to be due to 
the approximation that the photon number of the state |/3,C) is close to the photon number of the input state. The 
average value of this photon number is close to the photon number of the input state, however each individual value 
is not necessarily close to fi. The expression for the phase variance introduced depends on the inverse of the photon 
number, and the average of an inverse is not necessarily equal to the inverse of an average. The general expression is 



In figure [| we have also plotted the estimated theoretical limit based on the average of l/n p for the data obtained in 
the time dependent e case. Specifically, the expression plotted is 



As can be seen, the phase variance introduced converges to this far more closely than to the limit based on the photon 
number of the input state. 



We have verified that the same result for the increase in phase variance with time delay for the mark I case with 
simplified feedback is obtained by the full perturbation theory calculation as for the highly simplified calculation in 
Our numerical results also verify the extra phase variance of r/2 quite accurately. 

The result for the increase in the mark I phase variance only holds if the phase estimate at the end of the measurement 
is the final value of the intermediate phase estimate. If the actual value of arg A is used, the phase variance decreases 




(42) 




(43) 



6. CONCLUSIONS 
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for moderate time delays. This is because the worse intermediate phase estimate reduces the value of \B\, making 
arg A closer to the best phase estimate, arg C. 

We have also shown that the complete version of the simplified calculation in (J] to determine the increase in the 
variance of the mark II phase estimate does not give convergent results. An alternative technique shows that the 
theoretical limit to the introduced phase variance when there is a time delay of r is r/(8n). This is a factor of 4 
smaller than the result obtained in ||, though it is obtained for different limit conditions. 

The introduced phase variance converges to the theoretical limit in the three different cases with unsimplified feed- 
back considered. (These are the case with arg A v feedback, and the two cases with arg(A^C^ _e ) feedback introduced 
m §.) For the case with simplified feedback, however, the phase variance is far above the theoretical limit (around 
10 times). In this case, for large time delays the phase variance converges to the heterodyne phase variance, as the 
intermediate phase estimate becomes very poor. 

It is possible to correct the simplified feedback to reduce the phase variance greatly, but even this corrected feedback 
does not give results close to those for unsimplified a,rgA v feedback. This indicates that if there is any significant 
time delay in the system, it is better to use unsimplified feedback, even though the processing of the data is likely 
to introduce a larger time delay. This makes the improved arg(A^C^ _e ) feedback even more attractive, as one of the 
main reasons for using the simpler argA„ feedback was that it allows the simplified, analogue feedback circuit to be 
used. 
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